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Abstract 



We consider the Hamihonian dynamics of spherically symmetric Einstein 
gravity with a thin null-dust shell, under boundary conditions that fix the 
evolution of the spatial hypersurfaces at the two asymptotically flat infini- 
ties of a Kruskal-like manifold. The constraints are eliminated via a Kuchai^- 
type canonical transformation and Hamiltonian reduction. The reduced phase 
space r consists of two disconnected copies of R"^, each associated with one di- 
rection of the shell motion. The right-moving and left-moving test shell limits 
can be attached to the respective components of T as smooth boundaries with 
topology M^. Choosing the right-hand-side and left-hand-side masses as con- 
figuration variables provides a global canonical chart on each component of T, 
and renders the Hamiltonian simple, but encodes the shell dynamics in the 
momenta in a convoluted way. Choosing the shell curvature radius and the 
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"interior" mass as configuration variables renders the shell dynamics trans- 
parent in an arbitrarily specifiable stationary gauge "exterior" to the shell, 
but the resulting local canonical charts do not cover the three-dimensional 
subset of r that corresponds to a horizon-straddling shell. When the evolu- 
tion at the infinities is freed by introducing parametrization clocks, we find 
on the unreduced phase space a global canonical chart that completely decou- 
ples the physical degrees of freedom from the pure gauge degrees of freedom. 
Replacing one infinity by a flat interior leads to analogous results, but with 
the reduced phase space U M^. The utility of the results for quantization 
is discussed. 

Pacs: 04.20.Fy, 04.40.Nr, 04.60.Kz, 04.70.Dy 
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I. INTRODUCTION 



Spherically symmetric geometries have a long and useful history as a physically interest- 
ing and technically vastly simplified arena for gravitational physics. In vacuum, Einstein's 
theory with spherical symmetry has no local degrees of freedom, and the reduced phase space 
in the Hamiltonian formulation is finite dimensional. Including an idealized, infinitesimally 
thin matter shell brings in an additional finite number of degrees of freedom. Including a 
continuous matter distribution generically yields a (l+l)-dimensional field theory, with the 
exception of fields whose gauge symmetries exclude spherically symmetric local degrees of 
freedom. A familiar example of a field with such a gauge symmetry is the electromagnetic 
field. 

In this paper we consider spherically symmetric Einstein gravity coupled to an infinitesi- 
mally thin null-dust shell. From the spacetime point of view, the solutions to this system are 
well known (see, for example, Refs. fTj-Q), and they can be easily obtained from a junction 
condition formalism that is general enough to encompass null shells (see Ref. ^ and the 
references therein). Our purpose is to explore the Hamiltonian structure of this system, 
treating both the geometry and the shell as dynamical. Among the extensive previous work 
on Hamiltonian approaches to spherically symmetric geometries (for a selection in a variety 
of contexts, see Refs. fi-^), we follow most closely the canonical transformation techniques 



of Kuchaf |T0[ . Our main results can be concisely described as generalizing the spherically 
symmetric vacuum Hamiltonian analysis of Ref. [|T^] to accommodate a null-dust shell. 

Finding a suitable action principle requires care. The shell stress-energy tensor is a delta- 
distribution with support on the shell history, which is a hypersurface of codimension one. 
Einstein's equations for the system therefore admit a consistent distributional interpretation 



4^ , and the content of these equations is captured by the junction condition formalism of 
Barrabes and Israel 0. We recover these equations from a variational principle. We take 
the shell action to be that of a spherically symmetric thin cloud of radially-moving massless 
relativistic point particles, and we vary the total action independently with respect to the 
gravitational variables and the shell variables. We shall see that this variational principle 
can be made distributionally consistent and that the variational equations do reproduce 
the correct dynamics. Achieving this requires, however, a judicious choice of the regularity 
properties of the metric. 

We begin, in section by setting up the Hamiltonian formulation of the system in 
the Arnowitt-Deser-Misner (ADM) gravitational variables. The spacetime is taken to have 
Kruskal-like topology, with two asymptotically fiat infinities, and the spatial hypersurfaces 
are taken to be asymptotic to hypersurfaces of constant Killing time at each spacelike infinity. 
The Killing time evolution of the hypersurfaces is prescribed independently at each infinity. 
We specify the regularity properties of the gravitational variables, and demonstrate that the 
variational principle is consistent and leads to the correct equations of motion. 

In section |ITl| we perform a canonical transformation to a new chart in which the con- 
straints become exceedingly simple. Two of our new variables are Hawking's quasilocal mass 
M(r) and the two-sphere curvature radius R(r), just as in the vacuum analysis of Ref. []10 



However, to maintain a consistent distributional interpretation of the variables in the new 
chart, we are led to relate the momentum conjugate to M(r) to the Eddington-Finkelstein 
time whose constant value hypersurface coincides with the classical shell history, and not 
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to the Killing time as in Ref. |T0|. The momentum conjugate to R(r) needs to be modified 
accordingly. Remarkably, the canonical transformation can then be chosen to leave the shell 
canonical pair invariant. The transformation is mildly singular for geometries in which the 
shell straddles a horizon, but it can be extended to this special case in a suitable limiting 
sense. 

In section |^ we eliminate the constraints by Hamiltonian reduction. The reduced phase 
space r turns out to have dimension four. As the vacuum theory under our boundary 
conditions has a two-dimensional reduced phase space and as a test shell in a fixed 



spherically symmetric background has a two-dimensional phase space, this is exactly what 
one would have anticipated. We first obtain canonical coordinates (m_|_, m_, p+, p_) in 
which the configuration variables m± are the Schwarzschild masses on the two sides of the 
shell. The momenta p± can be interpreted as the Eddington-Finkelstein time differences 
between the shell and the infinities, after introducing an appropriate correspondence between 
our spatial hypersurfaces and hypersurfaces that are asymptotically null. The configuration 
variables m-|- are constants of motion, while the shell motion is indirectly encoded in the 
dynamics of p±. These coordinates become singular for horizon-straddling shells, but a 
global chart covering also this special case can be obtained by introducing suitable new 
momenta. We find that F consists of two disconnected copies of R^, each associated with 
one direction of the shell motion. The right-moving and left-moving test shell limits, in 
which the shell stress-energy tensor vanishes, can be attached to the respective components 
of r as smooth boundaries with topology M.^. 

In section ^ we introduce on F a local canonical chart in which the shell motion becomes 
more transparent. Assuming that the shell does not a straddle a horizon, the shell history 
divides the spacetime into the "interior," which contains a Killing horizon bifurcation two- 
sphere, and the "exterior," which does not. We choose the configuration variables in the 
new chart to be the curvature radius of the shell two-sphere and the interior mass, in an 
arbitrarily specifiable stationary exterior coordinate system. One can argue that this yields a 
Hamiltonian description of interest for an observer who scrutinizes the shell motion from the 
exterior asymptotic region, especially if the observer's ignorance of the interior asymptotic 
region is incorporated by setting the interior contribution to the Hamiltonian to zero. We 
give three examples of stationary exterior coordinate systems in which the Hamiltonian 
can be found in closed form. Also, choosing the spatially fiat exterior gauge [^3|-|i5[], and 



performing a partial reduction by setting the interior mass equal to a prescribed constant, we 
reproduce the spatially fiat shell Hamiltonian previously derived in Refs. pn|p7[ by different 
methods. 

In section |VT| we free the evolution of the spatial hypersurfaces at the spacelike infinities 
by introducing parametrization clocks. We find on the unreduced phase space a canonical 
chart in which the physical degrees of freedom and pure gauge degrees of freedom are com- 



pletely decoupled, in full analogy with the vacuum analysis of Ref. |Ty]. The pure gauge 
chart can be chosen so that the configuration variables are the curvature radius of the two- 
sphere and the Eddington-Finkelstein time, with the latter one appropriately interpreted 
across the horizons. 

In section |V11] we replace the Kruskal spatial topology S"^ x M by the spatial topology M^. 
The spacetime has then just one asymptotic region, and when the equations of motion hold, 
the spacetime interior to the shell is fiat. As in the Kruskal case, we take the asymptotic 
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region to be asymptotically flat, and we prescribe the evolution of the spatial hypersurfaces 
at the spacelike infinity. We then carry out the canonical transformation and Hamiltonian 
reduction. Expectedly, the reduced phase space turns out to consist of two disconnected 
copies of M^, with only the counterpart of the pair (m_|_, p_|_) of the Kruskal theory surviving. 
We conclude in section |V111| with a summary and a brief discussion, including remarks 



on the potential utility of the results in view of quantization. Some of the technical detail 
of the ADM dynamical analysis is postponed to the appendices. 

We work in Planck units, h = c = G = 1. Lowercase Latin tensor indices a,b, . . . 
are abstract spacetime indices. Dirac's delta-function is denoted by 6, while 6 denotes a 
variation. The curvature coordinates (T, R) for the Schwarzschild metric are coordinates in 
which the metric reads 

ds^ = -{l-2M/R)dT'^ + {l-2M/RY^dR^ + R^dn'^ , (1.1) 

where dVL^ is the metric on the unit two-sphere and M is the Schwarzschild mass. T and R 
are called respectively the Killing time and the curvature radius. 



II. METRIC FORMULATION 

In this section we set up the Hamiltonian formulation for spherically symmetric Einstein 
gravity coupled to a null-dust shell. We pay special attention to the regularity of the 
gravitational variables and the global boundary conditions. 

A. Bulk action 

Our spacetime geometry is given by the general spherically symmetric Arnowitt-Deser- 
Misner (ADM) metric 

ds^ = -N^dt^ + k^{dr + N^'dtf + R^dn^ , (2.1) 

where dVi^ is the metric on the unit two-sphere, and A^, A^*", A, and R are functions of the 
coordinates t and r only. Partial derivatives with respect to t and r are denoted respectively 
by overdot and prime, ' = d/dt and ' = d/dr. We take the spacetime metric to be 
nondegenerate, and A^, A, and R to be positive. 

The matter consists of an infinitesimally thin shell of dust with a fixed total rest mass m, 
which we initially take to be positive. Denoting the shell history by r = r(t), the Lagrangian 
matter action is 



Sl = -m dtdm-K^ix + m] , (2.2) 



where a hat is used to denote the value of a variable at the shell. The shell can be envisaged 
as a thin spherically symmetric cloud of radially-moving massive relativistic point particles. 

The Lagrangian gravitational action is obtained by performing integration over the angles 
in the Einstein-Hilbert action, (IGvr)"^ J d^x^J—g R. Discarding a boundary term, the result 
is [|,0Jig,|23,|g,|39| 
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dt / dr 



N-^ \r[A - (AA^'^)'](i? - R'N'') + iA(i? - R'N"" 



SI 



+N i^A-^RR'A' - A^^RR" - \A-^R!^ + 
The Lagrangian action of the coupled system is 

Sl = S^j^ + S\^ + boundary terms. 



(2.3) 



(2.4) 



We shall consider the regularity properties of the variables, the boundary conditions, and 
boundary terms after passing to the Hamiltonian formulation. 

The momenta conjugate to the configuration variables r. A, and R are 



N'i-A'^ix + N' 



Pa 
Pr 



-^{R-N^R') , 

-^iR-N''R')-§ A-{N^Ay 



A Legendre transformation gives the Hamiltonian bulk action [27,^] 
5e = J dt pi + J dr (^PaA + PrR - NH - A^"i7^) 



(2.5a) 

(2.5b) 

(2.5c) 

(2.6) 



where the super-Hamiltonian constraint H and the radial supermomentum constraint Hr 
contain both gravitational and matter contributions. In the limit m —>■ 0, these constraints 
take the form 



H 



API 



PaPr ^ RR" 



2i?2 R ' A 
Hr = PrR' - P;A - p6{r - V 



RR'A' 
A2 + 2A 



(2.7a) 
(2.7b) 



where rj := sign(p). From now on, we shall work exclusively in this zero rest mass limit, 
with the bulk action (|2.6|) and the constraints (|2.7] ) . As will be verified below, the shell then 
consists of null dust. 

The Hamiltonian constraint ( |2.7a ) is not different iable in p at p = 0. As we shall verify 
below, an initial data set with nonzero p cannot evolve into a set with p = We assume 

from now on that p is nonzero: this breaks the phase space into the two disconnected sectors 
?7 = ±1. The limits p — > 0± within each sector will be addressed in subsection p.V (J| . 



B. Local equations of motion 

In the presence of a smooth matter distribution, one can assume the spacetime metric 
to be smooth (C°°). In the idealized case of an infinitesimally thin shell, the metric can be 
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chosen continuous but not differentiable at the shell P, |42|^ , |47[] . The issue for us is to find 
smoothness assumptions that give a consistent variational principle. We wish to make both 
the action ( |2.6| ) and its local variations well defined and such that the resulting variational 
equations are equivalent to Einstein's equations with a null-dust shell. 

We follow the massive dust shell treatment of Ref. [^. In contrast to the case of a 
massive dust shell, we shall find that the smoothness conditions introduced in Ref 
make our null-dust variational principle fully consistent. 



37 



As in Ref. [0], we assume that the gravitational variables are smooth functions of r, with 
the exception that A^', {N^)', A', R', P/^, and Pr may have finite discontinuities at isolated 
values of r, and that the coordinate loci of the discontinuities may be smooth functions of t. 
All the terms under the r-integral in the action (|2.6|) are well defined in the distributional 
sense. The most singular contributions are the explicit matter delta-contributions in the 
constraints, and the implicit delta-functions in R" and P^. All these delta-functions are 
multiplied by continuous functions of r. The remaining terms are at worst discontinuous 
in r. The action is therefore well defined. 

Local independent variations of the action with respect to the gravitational and matter 
variables give the constraint equations 



H 

Hr 







and the dynamical equations 
'APa Pr 



A = N 
R = 



P2 R 



- — + (A^^A)' 



R 



Pa = - 



A^ 
2 



p2 



P' 

P2 I A 



A^'PP' 



A^ 
A 



P = P 



API_PaPr 

P3 P2 



R!_ 
X 



A^'P 
XT 



A 



(2.8a) 
(2.8b) 



(2.9a) 
(2.9b) 
(2.9c) 

(2.9d) 
(2.9e) 
(2.9f) 



As will be discussed in detail in appendix ^, our smoothness conditions imply that equa- 
tions (p.8|) and all save the last one of equations (|2.9|) have an unambiguous distributional 



interpretation. Equation (|2.91]) , on the other hand, is ambiguous: the right-hand side is a 
combination of spatial derivatives evaluated at the shell, but these derivatives may be dis- 
continuous. We need to examine the dynamical content of the well-defined equations, and 
the possibilities of interpreting equation ( p.91| ). 

A first observation from equation ( p.9e|) is that the shell history is tangent to the null 
vector 1°" whose components are 
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t = 1 , (2.10a) 
t = 7]Nk-^ ~N~^ . (2.10b) 

For 7] = 1 [rj = —1), l"" is the future null vector that points towards relatively larger 
(smaller) values of r. From the definition of the shell stress-energy tensor, dgSshea = 
1 1 y/~9 d'^xT''^ S{gab), we find 

The shell is therefore null, with positive surface energy but vanishing surface pressure 0. 
This confirms that the shell consists of null dust. 

All solutions to the spherically symmetric Einstein equations with a null-dust shell can 
be found from a sufficiently general junction condition formalism On each side of 

the shell, the spacetime is locally part of the extended Schwarzschild geometry. If the 
global structure of the spacetime is Kruskal-like, with two asymptotically flat infinities, 
there are only two qualitatively different cases. First, if the shell is not static, the junction 
is completely determined by continuity of the two-sphere radius at the shell. The motion is 
clearly geodesic in each of the two geometries, and the radius of the two-sphere serves as an 
affine parameter in either geometry. The spacetime is either that shown in figure |l|, or its time 
and/or space inverse. Second, if the shell is static, the junction is along a common horizon, 
and the masses must agree. The soldering is affine, meaning that the affine parameters 
along the horizon with respect to the two geometries are affinely related; however, as the 
stress-energy tensor of the shell is by assumption nonvanishing, the bifurcation two-spheres 
on the two sides do not coincide. The spacetime is either that shown in figure ^ or its time 
inverse. 

Now, away from the shell, equations ( |2.8| ) and ( |2.9| ) are well known to be equivalent to 
Einstein's equations. We shall investigate equations (|2.8|) and (|2.9|) at the shell in detail in 
the appendices. The result is that, when combined with the fact that the geometry is locally 
Schwarzschild on each side of the shell, the well-defined equations, (|2.8| ) and (|2.9a| )- (|2.9e|) . 



are equivalent to the correct null-dust junction conditions at the shell. They further imply 
that the right-hand side of (|2.9f|) is unambiguous, and that (|2.9fl) is satisfied as an identity. 



Our variational principle is therefore consistent, and it correctly reproduces the motion of a 
null-dust shell. 

A check on the consistency of our formalism is that the Poisson brackets of our constraints 



can be shown to obey the radial hypersurface deformation algebra [4S], as in the absence of 



the shell, and as with a massive dust shell We therefore have a Hamiltonian system 



with first class constraints |H9 . 



C. Falloff 



What remains are the global boundary conditions. We take the coordinate r to have the 
range (—00,00), and as r ^ ±00 we assume the falloff [p!0| , ^] 



A(t,r) = l + M±\r\'^ + 0' 



(2.12a) 
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R{t,r) 


= \r\ + 0°°{\r\-') 




= 0^{\rn , 


PrM 


= 0-^{\r\-'-') , 


N{t,r) 


= Ar^ + 0°°(|r|"' 


N''{t,r) 





(2.12b) 
(2.12c) 
(2.12d) 
(2.12e) 
(2.12f) 

where M± and N± are functions of t, and e is a parameter that can be chosen freely in the 
range < e < 1. Here, 0°° stands for a term that falls off as r ^ ±oo as its argument, and 
whose derivatives with respect to r and t fall off accordingly. These conditions imply that 
the asymptotic regions associated with r ±oo are asymptotically flat, with the constant t 
hypersurfaces asymptotic to hypersurfaces of constant Minkowski time. are the rates at 
which the asymptotic Minkowski times evolve with respect to the coordinate time t. When 
the equations of motion hold, M± are time-independent and equal to the Schwarzschild 
masses. 

In the variational principle, we take N± to be prescribed functions of t, but leave M± 
free. The appropriate total action then reads 



S = Sj: + S9j, , (2.13a) 

where the boundary action is 

S9s = - J dt {N+M+ + N_M_) . (2.13b) 

The global structure of the spacetime is Kruskal-like, with two asymptotically flat asymp- 
totic regions. The classical solutions under these boundary conditions are precisely those 
described above and shown in figures |I| and 0. 

III. CANONICAL TRANSFORMATION 

In this section we find a new canonical chart in which the constraints become exceedingly 
simple. Away from the shell, our treatment closely follows that given by Kuchaf in the 
llO| . The new elements arise mainly from patching the two vacuum regions 



vacuum case 



together at the shell. 

Our canonical transformation turns out to be mildly singular when the masses on the 
two sides of the shell agree and the shell straddles a common horizon. We first perform the 
transformation, in subsection |111 A|, assuming that this special case has been excluded. We 



then argue, in subsection |III B| , that the transformation can be extended to the special case 
in a suitable limiting sense. 



A. Shell not on a horizon 



In the vacuum theory, Kuchaf [|T0[ found a transformation from the canonical chart 



(A, R, Pa, Pr] to the new canonical chart (M, R, Pm, Pr) defined by 



9 



M 
R 

Pm 
Pr 



iRil-F) , 
R , 

Pr - IR-'APa - \R-^F-^APa 



-R-'A-'F-' [{APAyiRR') - {APa){RR')'] , 



(3.1a) 
(3.1b) 
(3.1c) 

(3.1d) 



where 



F :-- 



Pa 
R 



(3.2) 



When the equations of motion hold, M is independent of both r and t, and its value is just 
the Schwarzschild mass. Similarly, when the equations of motion hold, we have Pm = —T', 
where T is the Killing time. The vacuum constraints can be written as a linear combination 
of M' and Pr, and the dynamical content of the theory becomes transparent. 

In the presence of our null shell, the variables (M, R, Pm, Pr) become singular at the 
shell. To see this, consider a classical solution in which the shell history does not lie on a 
horizon. As M is discontinuous at the shell, M contains at the shell a delta-function in r. 
As Pm is discontinuous at the shell, the product PmM is ambiguous. One therefore does 
not expect the chart (M, R, Pm, Pr) to be viable in the presence of the shell. 

To overcome this difficulty, we keep M and R but replace the momenta by ones that are 
smoother across the shell. We define first [p!0[ 



and 



R^,Pa 
A R 



R' Pa 
— ± ri — 
A ' R 



(3.3) 



(3.4) 



Note that F = F^F_ = F^F_^. When the equations of motion hold, -F+ vanishes on the 
left going branch (es) of the horizon and F_ vanishes on the right going branch (es) of the 
horizon. It follows that F_^ is nonvanishing on the horizon that the shell crosses. Now let 



Pm + VF^'R! 

vAf:^] , 

Pr - r]F-^M' 



r]A 



Pn + r/P(ln |P_,|)' + ^ (P_, - Fl^) 



When the equations of motion hold, equation ( |3.5a|) shows that 



n 



M 



where r* is the tortoise coordinate 147 



-(T- 
For 1] 



rjr 



(3.5a) 
(3.5b) 

(3.6) 



— 1), T — rjr* is the retarded (ad- 
vanced) Eddington-Finkelstein time coordinate. While Pm was associated with the Killing 



+ 1 (r] 
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time |T0|, our prospective new momentum Um is therefore associated with the retarded or 



advanced Eddington-Finkelstein time. 

Away from the shell, a calculation of the Poisson brackets shows that the set 
{M,R,IIm,^r) is a candidate for a new gravitational canonical chart. We need to find 
shell variables that complete this set into a full canonical chart. 

For the rest of this subsection, we assume that the shell history does not lie on a horizon. 
The special case where the shell straddles a horizon will be discussed in subsection [111 lj| . 

As a preliminary, suppose that the constraints ( |2.8| ) hold, and consider the regularity of 
the variables. Away from the shell, the constraints ( |2.8| ) imply that R' and Pa are continuous, 
and F± are thus both continuous. In the notation of appendix the distributional content 



( |A^ ) of the constraints at the shell can be written as 



= AF_ 




P + A(APa) 



(3.7a) 
(3.7b) 



From ( |3.7a| ) we see that is continuous at the shell. Equation ( |3.5a| ) then implies that 11^ 
is continuous, with the exception that it diverges on the horizon that is parallel to the shell 
history. The first equality sign in ( |3.5b| ), and the observation that the vacuum constraints 
are linear combinations of M' and Pr 0], imply that IIr is vanishing everywhere except 



possibly at the shell. The rightmost expression in ( 3.5b| ) shows that IIr cannot contain a 
delta-function at the shell, and IIr is therefore everywhere vanishing. From now on, we can 
therefore proceed assuming that and IIr are continuous, and that their r-derivatives 
have at most finite discontinuities at isolated values of r. By (|3.5a| ), the same will then 
hold for IIm, with the exception of the horizon where Um diverges. This tightens the 
neighborhood of the classical solutions in which the fields can take values, but it will not 
affect the critical points of the action. The reason for this assumption is that it will make 
the terms UmM and IIrR in our new action distributionally well defined. 
We can now proceed to the Liouville forms. A direct computation yields 



PaSA + PrSR + MSUm - UrSR 



+d 



{r]RdR\n\F^.^\y 
>PA 



{FI^ - F^r,) +r]RR'\n\F_r,\ 



(3.8) 



The variation 5 affects the smoothness of the gravitational variables in the same way as the 
time derivative in subsection [11 B| . Away from the horizon parallel to the shell history, on 
which vanishes and Um diverges, all the terms in ( |3.8|) are therefore distributionally 
well defined: the terms on the left-hand side are at most discontinuous in r, while the terms 
on the right-hand side may contain at worst delta-functions arising from (dR)'. The status 
at the horizon on which vanishes will be discussed below. 

To obtain the difference in the prospective Liouville forms, we need to integrate the 
relation ( p. 81) over r. In an integral over a finite interval in r, the only subtlety arises from 



the horizon on which F_^ vanishes. On a classical solution with mass M, it can be shown 



from the embedding analysis of appendix that 



AM 



Th) + 0{{r - rhY) 



(3.9) 
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where the subscript h indicates the values of the quantities at the horizon on which 
vanishes. Equations (p.9|) and (p. 5a ) therefore show that, on the classical solution, the 



integral of ( |3.8| ) across the horizon is well defined in the principal value sense, just as in the 
corresponding analysis of Ref . [|l^ . To extend this argument off the classical solutions, we 



note that when the constraints hold, M(r) is constant in r across this horizon. As our action 
contains the constraints with their associated Lagrange multiplies, we argue that M(r) can 
be assumed smooth at the horizon in the relation ( |3.(j| ). We can then again employ ( p.9| ) 
and ( p.5a|) , and it is seen as above that the integral of (|3.8| ) across this horizon is well defined 
in the principal value sense. 

What needs more attention is the falloff in ( |3.8| ) at the infinities. From (p.l2|) , ( p.lD , 
and (|3.5| ), we have 

M{t,r) = M±{t) +0°°{\r\~') , (3.10a) 
R{t,r) = \r\ + 0^{\r\-') , (3.10b) 
UM{t,r) = ±r]{l + 2M±\r\-')+0^{\r\-'"') , (3.10c) 
nR(t,r) = 0°°(|rr^-^) . (3.10d) 

This means that the integrals of the third term on the left-hand side and the total variation 
term on the right-hand side diverge as r — *• ±cx). The geometrical reason for this divergence 
is, as seen from (|3.6| ), that Um is associated with a null time, rather than an asymptotically 
Minkowski time. 

The cure is to introduce convergence functions that provide the necessary translation 
between asymptotically spacelike hypersurfaces and asymptotically null hypersurfaces. To 
this end, let g{M^, M_;r) be a function that is smooth in r and depends on our variables 
only through M_|_ and M_ as indicated. Let g have the falloff 

g{M+,M_;r) = ±Ml\r\''^ + 0°^{\r\~^~') . (3.11) 



Adding —rjdg on both sides of ( |3.8[ ) yields now an equation whose both sides can be in- 
tegrated in r from — oo to oo. The substitution terms arising from the first term on the 
right-hand side vanish, and we obtain 

/oo /"OO 
{PaSA + PRdR)dr= / {UrSR - MSUm + V^g)dr 
-oo J — oo 

+ ^{/" ^{F~i-F-,)+vRR'i^\F-,\~V9 dr"^ . (3.12) 



All the terms in (|3.12| ) are well defined, provided the integral across the horizon on which 
vanishes is interpreted in the principal value sense. We therefore see that the set 
(M, R, r, XIa/, IIr, p) provides a new canonical chart on the phase space. Note that this 
canonical transformation leaves the shell variables (t, p) entirely invariant. The geometrical 
meaning of the convergence function g will be discussed in section [TV] . 

What remains is to write the constraint terms in the action in terms of the new variables. 
Consider first the constraints away from the shell. A straightforward rearrangement yields 

NH + N''Hr = N^Ur + NM' , (3.13) 
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where 



N :-- 



:= N^'R' - NR~^P^ 



(3.14a) 
(3.14b) 



Note that A^'^ is the same as in Ref. [0. Both terms on the right-hand side of ( |3.13|) are 
distributionally well defined. Away from the shell, we can therefore include the constraints in 
the action in the form shown on the right-hand side of ( p.l3|) , with and A^'^ as independent 
Lagrange multipliers. This constraint redefinition is mildly singular on the horizon parallel 
to the shell history, owing to the divergence of Hjv/; however, one can argue as in Ref. |1^ 
that the redefined constraints are equivalent to the old ones by continuity. The falloff of the 
new multipliers is 



A^ 



TAr± + 0°°(|r| 



(3.15a) 
(3.15b) 



To recover the delta-constraint ( ^.7aD , we observe from ( ^.5b| ) that ( |3.7a| ) is equivalent 
to IIr not having a delta-contribution at the shell. We therefore argue that including in the 
action the constraint term — J dt J^^drN^Uf^, with A^*^ an independent Lagrange multi- 
plier, yields both the constraint IIr 
at the shell .0 



away from the shell and the delta-constraint (|3.7a]) 



Finally, consider the delta-constraint ( p.7b| ). Using ( |3.5a] ) and ( p.7aD , equation ( |3.1a] ) 
implies 



AM 



A(APa 



n 



M 



When ( p.7a|) holds, (|3.7b|) is thus equivalent to 

P 



AM 



n 



(3.16) 



(3.17) 



M 



Including in the action the constraint term — J dt dr N[M' — pn^/(5(r — r)], with N an 
independent Lagrange multiplier, therefore yields both the constraint M' = away from the 
shell and the delta-constraint ( p.7b|) at the shell. 
These considerations have led us to the action 



S 



dt 



pi+ / dr[URR-MUM + m 



J dt J dr^N^UR + N[M'-pU]^^6{r-x)]^ 



dt {N+M+ + N_M. 



(3.18) 



subtlety in this argument is that A^'^ need not be continuous at the shell, not even on the 
classical solutions. The product A^'^IIr would therefore not be distributionally well defined in the 
event that IIr did contain a delta-contribution at the shell. 
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Both the action and its variations are well defined. The Poisson bracket algebra of the 
constraints clearly closes. Note that the convergence term rig in no way contributes to the 
local variations of the action. 

The Liouville term —MIIm can be brought to a form in which the time derivative is 
on M, at the cost of introducing another convergence term. Let G{r) be a smooth function 
of r only, with the falloff 

G{r) = ±l + 0^{\rf^'') . (3.19) 

We then have 

r]g- MIIm = iUM-vG)M -r]g + j^{T]GM + 27]g- MUm) • (3.20) 



All the terms in (|]20D are well defined, and each side can be integrated in r from — oo to oo. 



We arrive at the action 



S = j dt\pi + j dr [HrR + (Hm - vG)M - r?^] | 

- J dt J c/r I A^'^Hr + iV[M' -pn^/(5(r -r)]} 

- J dt {N+M+ + A^_M_) . (3.21) 
The geometrical meaning of the convergence terms will become explicit in section [TV] . 



B. Shell on a horizon 



In subsection [III A| we excluded the special case where the shell straddles a horizon. We 



now discuss how this special case can be included. 



When the shell straddles a horizon, the zero of occurs at the shell. The delta- 



constraints at the shell are given by ( p.7|) . When the equations of motion hold, the masses 
on the two sides agree, and the embedding analysis of appendix y shows that equation 
( p.9|) holds, now with A/^ = A and rh = t. From (|3.5b|) we see that IIr cannot contain a 
delta- function at the shell. We can therefore again assume that and IIr are continuous, 
and that their r-derivatives have at most finite discontinuities at isolated values of r. 

The new feature in equation ( |3.8| ) is that the singularity of FZ^ and IIm now occurs at 
the shell. When the equations of motion hold, we see from ( |3.5a|) and ( |3.9| ) that integrating 



each side of ( p.8|) in r across the shell is well defined in the principal value sense, and we argue 
as above that this conclusion can be extended away from the classical solutions provided 
the constraints are understood to hold. We argue similarly that the left-hand side of ( p.8|) 
contains no delta-contributions at the shell, and it is therefore justified to interpret the 
integral of (|3.^ ) over r as the principal value. Convergence at the infinities is accomplished 
as above, and the substitution terms from the total r-derivative on the right-hand side of 
]8D vanish. We therefore again arrive at ( p.l2|) . Equations (|3.16|) and ( p.l7| ) remain valid, 



with the understanding (^m^ = 0, and the delta-constraints can be taken in the action 
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as before. To justify the manipulations leading to the action ( p. 211) , we again appeal to the 



constraints to argue that M can be regarded as smooth in r at the shell, and that M then 
does not contain a delta-function at the shell. 

We therefore see that the actions (|3.18|) and (|3.21|) remain valid in a suitable limiting 
sense also for a horizon-straddling shell. 



IV. REDUCTION 

In this section we eliminate the constraints and find the dynamics in the reduced phase 
space. We shall continue to treat the cases i] = ±1 separately, and we denote the correspond- 
ing two components of the reduced phase space by F^. We first assume, in subsection |IV A| , 



that the shell history does not lie on a horizon, and we then include the horizon-straddling 
shell as a limiting case in subsection [IV B| . Finally, in subsection [IV Q , we attach the right- 



moving and left-moving test shell limits to the respective components of the reduced phase 
space as regular boundaries. 

It will be useful in the reduction to assume a more definite form for the convergence 
function g. From now on, we take 

g{M+,M_;r)=Mlg4r) + M'g4r) , (4.1a) 
where g±{r) are smooth functions of r only, with the falloff 

g±{r) = ±\r\-^e{±r) +0°°{\r\-^~') , (4.1b) 
where 9 denotes the step function. 

A. Shell not on a horizon 

In this subsection we assume that the shell history does not lie on a horizon. 
Solving the constraint IIr = implies that R and IIr simply drop out of the action. To 
solve the remaining constraint, 

UmM' - p6{r - r) = , (4.2) 

we write 

M = m+6'(r - r) + m_6'(r - r) , (4.3) 
where in±{t) are regarded as independent variables. We then have 

M(r) = m_|_^(r — t) + m_6'(r — r) + (m_ — m_|_)r5(r — r) , (4.4) 
and the constraint ( [4.2| ) implies 

p = (m+ - m_)n^ . (4.5) 



15 



Note that as the shell history does not lie on a horizon, each of the two factors on the 
right-hand side of ([4.5|) is nonvanishing. 



Using (|T51), (0), and (g^]), we find 



pr + 



dr 



(Hm - VG)M - r]g 



p+m+ + p_m_ 
d 



+ 



dt 



m 



Gdr 



where 



P^ 
P 



/oo 
dr [IlMd{r — t) — rjG6{r) — 2rjV[i^g^] 
-oo 

/oo 
dr [njvf6'(r — r) — rjG6{—r) — 2r]m_g_ 
-oo 



(4.6) 

(4.7a) 
(4.7b) 



The singularity of nM('") occurs in precisely one of the two integrals in (|4.71 ), and the integral 
over this singularity is interpreted in the principal value sense. Substituting ( [4.6| ) into ( p. 211 ), 
and dropping the integral of a total derivative, we obtain the reduced action 



S = J dt (p+m+ + p_m_ — A^_|_m_|_ — A^_m_) . 



(4.8) 



This shows that the set (m_(_, m_, p+, p_) provides local canonical coordinates on f^. The 
equations of motion derived from the action ([4. 8]) read 



P± 



, 



(4.9a) 
(4.9b) 



The emergence of m-|- as two coordinates on is not surprising: on a classical solu- 
tion, m± are the two Schwarzschild masses, and these masses together with r] completely 
determine the four- dimensional spacetime. To understand the geometrical meaning of p±, 
we recall from (|3.6| ) that without the convergence terms proportional to G and g±, the 
integrals in ( |4.7| ) would give the Eddington-Finkelstein time differences between the shell 
and the infinities on the constant t hypersurface. As the constant t hypersurface extends 
to the spacelike infinities, such null-time differences would be infinite. The role of the con- 
vergence terms in ( [4.7|) is to absorb the infinities: one can think of the convergence terms 
as associating to the constant t hypersurface a hypersurface that is asymptotically null as 
r ±oo. For rj = —1, this associated hypersurface extends from the left-hand-side T"*" 
to the right-hand-side I~ (points A and B in figure ||); for r] = 1, the situation is the 
reverse. Thus, p_|_ is the Eddington-Finkelstein time difference between the shell and the 
right-hand-side infinity of the associated asymptotically null hypersurface, and — p_ is the 
Eddington-Finkelstein time difference between the shell and the left-hand-side infinity of the 
associated asymptotically null hypersurface. The equations of motion (|4.9b|) show that the 
time evolution of p± only arises from the evolution of the constant t hypersurfaces at the 
infinities. Thus 
is encoded in equations ( [4.9 b| ) 



in this local canonical chart on F^, the information about the shell motion 
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It should be emphasized that the degrees of freedom present in p± are invariant under the 
isometries of the spacetime. Kilhng time translations on the spacetime move both the shell 
history and the constant t hypersurface: in particular, they move the two asymptotic ends of 
the constant t hypersurface, and hence the asymptotic ends of the associated asymptotically 
null hypersurface. However, the Eddington-Finkelstein time differences that constitute the 
momenta are invariant under Killing time translations. 

As we have assumed that the shell history does not lie on a horizon, the coordinates 
(m+, m_, p+, p_) do not form a global chart on F.^. Instead, these coordinates provide two 
disjoint local canonical charts, covering two disconnected sets in F^: one for < m_ < m_|_ 
and the other for < m_|_ < m_, with unrestricted values of p± in each chart. These 
coordinates cannot be extended to m+ = m_. The reason is that when m_|_ = m_, the 
shell history lies on a horizon, the singularity in IIm is at r = r, and the first term under 
each integral in ( [4.7|) makes both p_|_ and p_ divergent. We shall address the special case 
m_|_ = m_ and the global structure of F^ in subsection IV B. 



It will be useful to introduce on F^ another set of local canonical coordinates, 
(m, m, p, p) , by the transformation 

m=|(m+ + m_) , (4.10a) 

m = |(m+ - m_) , (4.10b) 

P = P+ + P- , (4.10c) 

p = p+-p_ . (4.10d) 



The inverse transformation is 



m-t = m ± m , (4.11a) 
p± = i(p±p) . (4.11b) 



From (^4.7|) , (|4.10c| ), and ( 4.10d ), we see that p contains the information about the asymptotic 
ends of the constant t hypersurface, whereas the information about the location of the shell 
with respect to the infinities is encoded in p. We can therefore loosely regard the pair (m, p) 
as describing the vacuum spacetime dynamics, and the pair (m, p) as describing the shell. 
While not literal, this view will be helpful for understanding the global structure of F^ in 



subsections IV B| and IV C 



As defined by the transformation (|4.10|), the coordinates (m, m, p, p) provide two disjoint 



local canonical charts that cover on F^ the same two disconnected sets as the coordinates 
(m+, m_, p_|_, p_). The ranges of the variables in these two charts are respectively < m < 
m and < — m < m, each with unrestricted p and p. The coordinates (m, m, p, p) cannot, 
however, be extended to m = 0. While p remains finite for a horizon-straddling shell, it is 
seen from ( [4.7|) that p must diverge. 

B. Shell on a horizon 

We now wish to find on F^ coordinates that extend to the horizon-straddling shell. We 
shall first rely on the spacetime picture to identify the geometrical information that the 



17 



coordinates must carry in this limit. We then construct on a global canonical chart that 
contains this information. 

Consider the spacetime of figure [1|. The shell is left-moving, corresponding to 77 = —1, and 
the shell history lies in the future of the left-going horizon, corresponding to m+ > m_. The 
points A and B indicate the ends of the asymptotically null hypersurface that is associated 
to the hypersurface of constant t. is the difference in the Eddington-Finkelstein time 
between points pi and qi, and p_ is the difference in the Eddington-Finkelstein time between 
points q2 and p2. 

In this spacetime, let 71 be the radial null geodesic connecting pi to qi, let 72 the radial 
null geodesic connecting q2 to p2, and let 73 be the radial null geodesic connecting pi to p2- 
Let Aj {i = 1,2, 3) be the affine parameters on these geodesies, each normalized to have the 
range [0, 1]. A2 and A3 increase toward the future. Ai increases toward the future if qi is in 
the future of pi as shown in the figure, corresponding to p+ < 0, and it increases towards 
the past if qi is in the past of pi, corresponding to p+ > 0. In the special case p+ = 0, 
qi and pi coincide, and 71 degenerates to a point. We now define the quantities Q± by 



:= {d/dX,Ud/dX,)\ 
g- := id/dX2)Sd/d\s)X^ 



(4.12a) 
(4.12b) 



Similarly, consider a spacetime in which 1] = —1 but the shell history lies in the past of 
the left-going horizon, corresponding to m_|_ < m_. In this spacetime, the counterparts of 
points pi and p2 are below the left-going horizon, but the three null geodesies 7^ {i = 1,2,3) 
can be defined as above, the only modification being that the potentially degenerate one is 
now 72. In this spacetime, we again define Q± by ( [4.12|) . 

A straightforward calculation yields 



Q- 



5mj 



5 m 



m 



mi 



m+ + I m+ — m_ | exp 



P+ 



4m I 



m_ + mi 



m_ I exp 



P 



4m 



(4.13a) 
(4.13b) 



valid both for m+ > m_ and m_|_ < m_. The canonical coordinates (m+, m_, p+, p_) can 
therefore be replaced by the noncanonical coordinates (m+, m_, Q-), with < for 
m_|_ > m_ and Q+ < for m_|_ < m_. 

The crucial observation is now that Q±, as defined in (|4.12|) , remain well defined also for 
the spacetimes shown in figure ^ in which m_|_ = m_ and the shell history lies on a common 
horizon. In these spacetimes, Q± can each take arbitrary negative values. As the soldering 
along the horizon is affine, Q± precisely encode the coordinate-invariant information about 
the relative loci of the points pi, p2. A, and B (or, equivalently, the points pi, p2, qi, and 52)- 
This means that the set (m+, m_, Q^, QJ) provides a global, noncanonical chart on F_. The 
domain is Q_ < for m+ > m_, < for m_|_ < m„, and Q± < for m_|_ = m_. 

In the above construction we have taken r] = —1. It is clear that an entirely analogous 
discussion carries through for r] = 1, with straightforward changes in formulas ( |4.13| ), and 
yielding a global, noncanonical chart on on F+. 

To find a global canonical chart on F^, consider the transformation 
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TT := p + 877m (In |m/m| — 1) 
TT := p + Sr/m (In |m/m| + 1) 



(4.14a) 
(4.14b) 



Equations ( 4.14 ) clearly define a canonical transformation from (m, m, p, p) to (m, m, tt, tt) 
individually in the domains < m < m and < — m < m. It is straightforward to verify 



that the chart (m, m, tt, tt) 
values of tt and tt. For rj = 



becomes global on r„ when extended to m = with unrestricted 



■1, in particular, ( |4.13|) shows that Q± can be written as 



= — 16(m + m) < mexp 



m — m 



m + m 



TT + TT 



111 + m 



= — 16(m — m) < mexp 



m + m\ TT — TT 

m — m y 8(m — m) 



m 



m 



(4.15a) 
(4.15b) 



from which the regularity of the m — > limit is manifest. 

We have thus shown that the set (m, m, W, tt) provides a global canonical chart on F^. 
The domain of the variables is |m| < m, with W and tt taking all real values. We therefore 
^. The Hamiltonian reads 



have F^ ~ 



h = (A^+ + N_)m + {N+ - N_)m 



(4.16) 



The values of rn and m are constants of motion, whereas the equations of motion for W and 
TT show that the evolution of Q± only arises from the evolution of the constant t hypersurface 
at the two spacelike infinities. This means that the information about the shell dynamics 
is contained in the momentum equations of motion both for m 7^ (as was already seen in 
subsection |1V A|) as well as for m = 0. 



C. Test shell limit 

We have so far assumed that the unreduced shell momentum p is nonvanishing. We saw 
that this assumption is compatible with the dynamics, and that it divides the reduced phase 
space into the two disconnected sectors F^, labeled by r] = sign(p). As the unreduced bulk 
action (|2.6| ) is not differentiable in p at p = 0, it is not clear whether F_|_ and F_ are joinable 
to each other in any smooth sense. Our reduction formalism is not well suited to examining 



this issue: the canonical transformation of section |T| was tailored to the null hypersurfaces 
separately for 77 = ±1. 

We can, however, address the limit p — > individually in F+ and F_. As the shell 



stress-energy tensor (|2.11|) vanishes for p ^ 0, this is the limit of a test shell that traverses 
the spacetime without affecting it gravitationally. We shall now show that one can attach 
the right-moving and left-moving test shell limits respectively to F+ and F_ as smooth 
boundaries with topology M'^. 

When the test shell history does not lie on a horizon, the situation is straightforward. 
We can start with the coordinates (m_|_, m_, p_,_, p_), separately for < m_ < 111+ and 
< m_|_ < m_, and simply take the limit m+ = m_ with p± remaining finite. From 
the geometrical interpretation of p-t it is seen that this attaches to F„ those test shell 
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configurations in which the test shell does not straddle a horizon. The locus of the test shell 



history is determined by p± exactly as in subsection IV A . 

Including a horizon-straddling test shell is more intricate. In the global chart 
(m, m, 7f, tt), the limit of a test shell on a horizon is achieved by setting first m = and 
then taking 7]tt — oo while keeping W finite. On the other hand, the limit of a test shell off 
the horizon requires taking simultaneously m ^ and rjiT ^ — oo so that p and p remain 
finite. What we need is a new canonical chart in which both of these limits are brought to 
finite values of the coordinates. 

To this end, let 

x:=exp(^^^ , (4.17a) 

— = ! (4.17b) 
8m / 

n:=7f-™ . (4.17c) 
m 

As 

Tldm + PxSx = WSm + ndih — S{nm.) , (4-18) 



equations ([4.17|) define a canonical transformation from the chart (m, m, tt, tt) to the chart 



(m, X, n,pa;). The new canonical chart is global: the range is m > and x > 0, with 
unrestricted H and px- The qualitative location of the shell history is governed by the sign 
of Px'- Px > {px < 0) yields a shell in the future (past) of the horizon that is parallel to the 
shell history, while Px = yields a shell history on the horizon. It is now easily seen that in 
this chart the test shell limit is x ^ 0, with the other coordinates remaining at finite values. 
A horizon-straddling test shell is recovered with px = 0, whereas px > {px < 0) gives a 
test shell in the future (past) of the horizon that the test shell does not cross. Clearly, the 
test shell limit constitutes a smooth boundary of with topology M^. 



V. HAMILTONIAN FOR THE SHELL RADIUS IN STATIONARY EXTERIOR 

COORDINATES 



While the charts on introduced in section |I^ are well adapted to the geometry of the 
spacetime, they contain the information about the shell motion in a nontransparent manner. 
In this section we introduce on F^ a local canonical chart that describes more directly the 
motion of the shell in the spacetime geometry. A chart of this kind is of particular physical 
interest if one wishes to quantize the system as a model of black hole radiation with back 
reaction |0|,|l],||]. 

The physical situation we have in mind is a static observer who scrutinizes the shell 
motion from an asymptotically fiat infinity. For definiteness, we take this infinity to be 
the right-hand-side one. We set = 1, so that the coordinate time t coincides with the 
observer's proper time. To incorporate the observer's ignorance of what is happening at the 
left-hand- side infinity, we set = 0. 
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We further assume that the shell history reaches a future or past null infinity on the 
right-hand side. The Penrose diagram for 77 = —1 is therefore as in figure |l|, and the Penrose 
diagram for 77 = 1 is the time inverse. In particular, we have m_ < m+, and we are in the 
region of P^ covered by the chart (m+, m_, p+, p_) with < m_ < m+. 

Consider thus the chart (m-|., m_, p+, p_) with < m_ < m+.From section ^ we recall 
that the pair (m+,p+) only carries information about the geometry right of the shell, and 
the pair (m_, p_) only carries information about the geometry left of the shell. To describe 
the motion of the shell as seen from the right-hand-side infinity, we can therefore leave the 
pair (m_,p_) intact and seek a canonical transformation that replaces (m_|_,p_|_) by a new 
pair. 

To specify the new pair, we choose in the Kruskal geometry right of the shell a stationary 
coordinate system that conforms to the falloff ( |2.12| ) with = 1. As R{r) is then an 
increasing function, we can assume R{r) = r without loss of generality. Let f stand for the 
shell curvature radius in these coordinates: f(t) := R(t(t)) = t(t). We now seek a momentum 
p such that there is a canonical transformation from the pair (m_(_,p_|_) to the pair {f,p). 
Writing the transformation as p+ = p+(r, m^) and p = p{f, m_|.), the canonicality criterion 
reads 

d[p{r,m+)] ^ g[p+(r, m+)] 
9m_|_ df 

Substituting p+ from ( [4.7aD to the right-hand side of (lOP yields 

d[p{r, m+)] 



dm. 



nA/(r,m+) , (5.2) 



where Hiviir, m+) is determined by the choice of the stationary coordinate system. Note that 
the convergence functions G{r) and g+{r) have not entered (|5.2|) . Solving the differential 
equation ( |5.2|) for p{f, m+) yields the desired canonical transformation, and inverting this 
solution gives m+(f,p) as a function in the new canonical chart. The action reads 

S= (it [p_m_ + — m+(f,p)] . (5.3) 



The shell stress-energy tensor in the new chart can be found using ( |2.10| ), ( |2.11| ) and ( ^3| 



As explicit examples, we now present the shell Hamiltonians m+(f,p) in four different 
stationary coordinate systems. We arrived at the first three coordinate systems by seeking 
a simple functional form for m_|_(f,p). The fourth coordinate system is the spatially fiat one 
used in Ref. IBDI . 



A. Polynomial gauge 

As a first example, we consider coordinates in which the metric reads 

R = r , (5.4a) 

A^ = ^ ^ ^ 2m+/r + {2m+/rf + (2m+/r)^ , (5.4b) 

A^AT^ = -r]{2in+/rf . (5.4c) 
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With r > 0, these coordinates cover half of Kruskal manifold in the appropriate manner, 
and the falloff (|2.12 ) is satisfied with e = 1. The relation to the curvature coordinates is 



R = r , (5.5a) 

T = t + 2?7m_|_ In |1 — 2m_|_/r| . (5.5b) 

We find 

nM(r, m+) = 7^(1 + 2m+/r) . (5.6) 

Solving ( |5.2|) with a convenient choice for the integration constant, we obtain 

m+(r,p) = ^/r]pf — . (5.7) 

The equation of motion derived from the Hamiltonian ( |5.7| ) can be integrated as rjt = 
r + 2m+ ln(f/m_|_) + constant. It is easily verified that this is the correct equation for a null 
geodesic in the metric ( ^.4|) . 

B. Exponential gauge 

Consider next coordinates in which the metric reads 

R = r , (5.8a) 
= exp(2m4./r)[2 — (1 — 2m+/r) exp(2m_|_/r)] , (5.8b) 
A^N'' = -7^ [1 - (1 - 2m+/r) exp(2m+/r)] . (5.8c) 

With r > 0, these coordinates cover half of Kruskal manifold in the appropriate manner, 
and the falloff (|2.12|) is satisfied with e = 1. The relation to the curvature coordinates is 



R = r , (5.9a) 

/r 

(1 - 2m+/r'y^ - exp(2m+/r') dr' . (5.9b) 

We find 

IlAf (r, m+) = exp(2m+/r) . (5.10) 
Solving ( [5. 21 ) with a convenient choice for the integration constant, we obtain 

m_(_(f,p) = |f ln(2?7p/f) . (5-11) 
The equation of motion can be solved implicitly in terms of the exponential integral function. 
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C. Eddington-Finkelstein type gauge 



Consider next coordinates in which the metric reads 

R = r , (5.12a) 
A2 = = 1 + 2m+/r , (5.12b) 
K^N" = -2r]m+/r . (5.12c) 

With r > 0, these coordinates cover half of Kruskal manifold in the appropriate manner. 
The relation to the curvature coordinates is 

R = r , (5.13a) 
T = t + 2r7m+ln|r/(2m+) - 1| . (5.13b) 

We recognize these coordinates as simply related to the Eddington-Finkelstein coordinates 



47| : t — rjr is the retarded (advanced) Eddington-Finkelstein time for rj = 1 (rj = —1). In 
terms of the tortoise coordinate r* := r + 2m+ ln[r/(2m_|_) — 1], we have t — rjr = T — rjr* . 

There is a minor technical issue in that the coordinates ( ^.12| ) do not obey the 
fallofF (|2.12| ): we have Pa = — 2r7m+(l + 2m+/r)^^^^, which violates (|2.12c|) . We there- 



fore take the coordinates ( [5. 12] ) to hold for r < -Rcut? where i?cut is a large parameter, and 
smoothly deform them to a faster falloff for r > i?cut- As equation ( |5.2| ) is local in f, the 
form of the canonical transformation for r < i?cut is independent of that for r > -Rcut- In 
the end, we can either leave the Hamiltonian unspecified for r > -Rcut, or argue that one can 
take the limit i?cut — > cxd in the sense of some suitable renormalization in the parameters of 
the canonical transformation. 

Proceeding in this way, we find 

nM(r,m+) =r/ , (5.14) 

and, with a convenient choice of the integration constant, 

m+(f,p)=r7p . (5.15) 

The Hamiltonian ( |5.15|) clearly correctly reproduces the fact that the Eddington-Finkelstein 
time t — rjr is constant on the shell history. 

D. Spatially flat gauge 

As the last example, we consider coordinates in which the metric reads 

R = r , (5.16a) 
A = A^ = 1 , (5.16b) 
A^'' = -r]^2ui+/r . (5.16c) 

With r > 0, these coordinates cover half of Kruskal manifold in the appropriate manner. 
The relation to the curvature coordinates is 
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R 



T = t + 2rj \ A/2m+r + m_|_ In 



(5.17a) 

(5.17b) 



1 + ^2^ 

We recognize these coordinates as the spatially flat coordinates [^3|-|45|], recently employed 



in the study of Hawking radiation with back reaction in Ref. ||30 



There is again a minor technical issue in that the coordinates ( p. 16] ) do not obey the 
falloff ( |2.12|) . A Hamiltonian falloff analysis compatible with these coordinates has been 



discussed in the metric variables in Ref. jS^]. Here, however, we shall simply argue in terms 
of a cutoff parameter Rent as above. We have 



HM(r, m_, 



1 + V2m+/r 



(5.18) 



and solving (|5.2| ) with a suitable integration constant yields m+(f,p) implicitly as the solu- 
tion to 



Tjp = A/2m_|_r — f In ^1 + \/2m+/r^ . 
In order to make a connection to the work in Ref. ||3^, we define 



Pc := P - 



V 



2m_f — f In ^1 + A/2m_/ 
- 2A/2m_f + 4m_ In (l + A/f/(2m_) 



(5.19) 

(5.20a) 
(5.20b) 



Equations ( ^.20| ) define a canonical transformation from the chart (f, m_,p, p_) to the new 
canonical chart (f, m_,pc, Pc). The Hamiltonian m+(f, m_,pc) in the new chart is obtained 
(in implicit form) by eliminating p from (|5.19|) and (|5.20a| ). As the value of m_ is a constant 
of motion, the system can be partially reduced by regarding m_ as a prescribed constant. 
The term J dt Pcrii- then drops from the action, and we obtain 



S = j dt[pcf - iLn+{f,m_,Pc)] . 



(5.21) 



For 1] = —1, this is the action derived in Refs. by different methods. For r/ = 1, it 

is not. The reason is that the coordinates ( |5.16 ) are the ingoing spatially fiat coordinates 
for 1] = —1 and the outgoing spatially fiat coordinates for rj = 1, thus covering all of the 
spacetime right of the shell in each case, whereas Ref. |]5D| was physically motivated to use 
the ingoing spatially fiat coordinates irrespectively the direction of the shell motion. It 
would be straightforward to repeat the above analysis with the sign of rj in (|5.16|) reversed, 
recovering the result of Refs. [3C,37] for t] = 1. Note, however, that with the sign of t] in 
( ^.16| ) reversed, the coordinates do not cover the part of the shell history that lies inside the 
horizon. 
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VI. PARAMETRIZATION CLOCKS AT THE INFINITIES 



In the previous sections we fixed the evolution of the spatial hypersurfaces at the spacehke 
infinities by taking N± to be prescribed functions of t. In this section we free this evolution 
by making the replacement [l^ 



A^4 



±f4 



(6.1) 



in the boundary term in the actions ( p. 13 ) and ( p. 21 ). The variations of N become then 
unrestricted at r ±oo, but varying the action with respect to t± yields the relations 
( |6.1| ) as equations of motion. The new variables t± are the proper times measured by static 
standard clocks at the respective infinities, with the convention that r+ increases toward the 
future and r_ increases toward the past. 



In the absence of a shell, it was shown in Ref. [10| that the action containing t± as 
independent variables can be brought to a canonical form in which the unconstrained degrees 
of freedom and the pure gauge degrees of freedom are entirely decoupled. We now outline the 
analogous result in the presence of the null shell. For brevity, we shall refrain from explicitly 
spelling out the smoothness properties of the various emerging phase space functions. 

We start from the action ( 3.21 ), with g given by ( [4.1|) , and we make in the boundary term 
the replacement ( p.l] ). The resulting action is a sum of two decoupled parts: a Hamiltonian 
action Sr consisting of the terms that contain the pair (R,nR), and the remainder 5*0. We 
only need to consider 5*0. As in section |ITT|, we assume first that the shell history does not 
lie on a horizon, and relax this assumption at the end. 

Under the time integral in 5*0, the terms homogeneous in the time derivatives are 



e := pir - M+f+ + M_f_ + 



dr 



[UM-vG)M-r]g 



(6.2) 



We pass from the noncanonical chart (r, M(r), p, nA^(r); r_|_, r_) to the new chart 
(m, r, r(r),p, p, nr(r)), defined by 



r(r) := M'(r) , 

nr(r) := i(r+ + r_) + [ dr' h[UM{r') - r]G{r')] [e{r' - r) - e{r - r')] 

J —oo 

-T]M+g+{r')+r]M^g4r')^ , 

m:=i(M+ + M_) , 

p := r+ - r_ + / dr UuM{r) - r]G{r)] - 2r]M+g+{r) - 2r]M.g^{r)\ 

J — oo 

The falloff is 

r(t,r) = 0-(|r|-^-^) , 
nr(t,r) = -2r/M±ln|r/M±| +0°°(|rn . 



(6.3a) 

(6.3b) 

(6.3c) 

(6.3d) 



(6.4a) 
(6.4b) 



By techniques similar to those in Ref. 



we find 
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e = pr + pm+ / rfrnr(r)r(r) + ^ (M+r+ + M_r_ 



(6.5) 



The chart (m,t,T{r),p,p,Ilr{r)) is therefore canonicaL Dropping the integral of a total 
derivative, the action reads 



5*0 



dt \pi + pm+ / drUrt 



drN 



- / dt 



(6.6) 



To express the constraint in the new chart, we have used equation ( |6.3aD and the relation 
Um = TjG — n'p, which follows by differentiating ( |6.3b| ). 

The action (|6.6| ) is canonical, but the constraint couples the variables in a nontransparent 
way. To decouple the degrees of freedom, we pass to the chart (m, m, r(r), p, p, np(r)), 
defined by 



r(r) := r(r) - p[r]G{r) - n^(r)] ^6{r - x 



Hf (r) := nr(r) 



-1 



m := Ip(7]G-U'^ 
p -=21^^-27] I drG{r) . 



The falloff of F and lip is clearly the same as that of F and Fir, given in 
analogue of relation ( |X^ ) for flp, we find 



pi-pm+ I dr [ HpF - HpF ) = ^ 



2r]m / drG{r) 



(6.7a) 
(6.7b) 

(6.7c) 
(6.7d) 

Using the 
(6.8) 



The chart (m, m, F(r), p, p, np(r)) is therefore canonical. Dropping the integral of a total 
derivative, the action reads 



So = dt [ pm + pm + / dr FIpF 



dt 



drNT . 



(6.9) 



The unconstrained canonical degrees of freedom, {m,m,p,p), have now become decoupled 
from the pure gauge degrees of freedom. 

To put the action in a more transparent form, we write 



V{r) ■=Uf,{r) + lp-ri [ dr' G{r') 

Jo 

dr' [UMir')e{r' - r) - 7]G{r')e{r') - 2T]M+g+{r')] 



= r+ + 
Uvir) := -f(r) 



(6.10a) 
(6.10b) 
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The falloff is 



V{t,r) = -r]\r\~2r]M±\n\r/M±\+0°^{\rf) , (6.11a) 
Uv{t,r) = 0°°{\r\-^-') . (6.11b) 



As 



/ dr (Uf.T -UvV) = -^p drliy + — drlifV, (6.12) 

the transformation to the chart (m, m, V"(r),p, ny(r)) is not canonical as it stands. How- 
ever, it becomes canonical after the first term on the right-hand side of (|6.12| ) is absorbed 
into the constraint term by writing 

N^-=-N + \f (6.13) 

and regarding A^^ as a new Lagrange multiplier. As the equations of motion imply p = 0, 
the falloff of A^^ is 

= ±N± + 0'^{\r\~') . (6.14) 
Dropping the integral of a total derivative, and including S'r, we finally obtain the action 

+ J dt J dr (UvV + UfiR- N^Uv - N^Ur^ . (6.15) 

All the variables in the action (|6.15|) have a transparent geometrical meaning. 
From (|6.3a|) , (|6.3c| ), (|6.7a|) , and ( |6.7c|) , we see that m and m are respectively equal to 



the variables m and fn introduced in section Similarly, using (|6.1| ), we see that p and 
p can be interpreted as the time-independent initial values of the variables p and p intro- 
duced in section As for the pure gauge degrees of freedom, R is the curvature radius, 
and equation (|6.10a]) shows that V is the Eddington-Finkelstein time. The action (|6.15 



therefore provides a natural generalization of the vacuum action given in Eq. (149) of Ref. 

We have here assumed that the shell history does not lie on a horizon. This assumption 
can be relaxed, in a suitable limiting sense, by performing on the coordinates (m, fh, p, p) 
transformations analogous to those given for the coordinates (m, m, p, p) in section 



VII. M? SPATIAL TOPOLOGY 

In this section we consider the canonical transformation and Hamiltonian reduction for 
spatial topology M^. For concreteness, we take the evolution of the spatial hypersurfaces at 
the single spacelike infinity to be prescribed as in section |T[ It will be seen that the reduced 
phase space consists of two disconnected components, one for an expanding shell and the 
other for a collapsing shell. Each component has the topology M^. 
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We start from the action principle. In the bulk action ( |2.6| ), we take < r < oo, with 
oo. As r 



the falloff (Pl^) as r 



A{t,r 
R{t,r 
PA{t,r 
PRit,r 
N{t,r 



0, we introduce the falloff 

= Ao + 0(r2) , 
= Rir + 0{r^) , 

= P^^r + 0{r^) , 
= N, + 0{r^) , 
= N{r + 0{r^) , 



(7.1a) 
(7.1b) 
(7.1c) 
(7.1d) 
(7.1e) 
(7.1f) 



where Aq > 0, -Ri > 0, Pp^.^, Pr^, Aq > 0, and A'[ are functions of t only. It is straightforward 



to verify that the falloff ( [7.1|) is consistent with the constraints and preserved by the time 
evolution. By ( 3.1a| ) and ( |3.2| ), the falloff ( [7.1| ) implies that the mass left of the shell must 
vanish when the equations of motion hold: r = is then just the coordinate singularity 
at the center of hyperspherical coordinates in flat space. The classical solutions therefore 
describe a shell with a flat interior, and the spatial topology is . The action appropriate 
for fixing A^_|_ is 



S = Sy 



with < r < oo. 



dt N,M. 



(7.2) 



where St. is given by 

The canonical transformation of section |T| goes through with the obvious changes. The 
new action is as in ( |3.18D , except that the integral is from r = to r = oo and the term 
N_M_ is missing. G{r) and g{M^;r) are smooth in r and have the same behavior as 
r +00 as in section |ITl| . The falloff of the new fields as r — can be found from ( [7.1D ; 
for example, we have 



N{t,r) = -NoAoRi' + 0{r) 



(7.3) 



All the new fields remain regular as r ^ 0. In particular, M{r) tends to zero as r ^ 0. 

The Hamiltonian reduction proceeds as in section with the simplification that the 
interior mass vanishes. The reduced phase space consists again of two disconnected compo- 
nents, denoted now by T^. We take g{M+; r) to be as in (|4.1|) with M_ = 0, and we solve 
the constraint ( [4. 2] ) as in ( [4. 3] ) and ( [4. 51 ) with m_ = 0. Note that as 11^/ has the same sign 
as p, equation ([4.5| ) implies m_|_ > 0. We find 



fOO 

pr + / dr 





(Um - r]G)M - r]g = p+m+ + ^ ( 77111+ / G dr 



where 



POO 

p+ := / dr [IlM9{r - t) - r]G - 2r]m+g+] 
Jo 



(7.4) 



(7.5) 



Substituting this in the action and dropping the integral of a total derivative, we obtain the 
reduced action 
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S = dt (p_|_m+ — A^+m+) 



(7.6) 



Thus, the pair (m_|_, p_,_) provides a canonical chart on T^. As Um does not have singularities, 
the definition ( [7. 51) is always good: the chart is global, and the topology of is M^. The 
test shell limit can be attached as a smooth boundary with topology M at m_|_ = 0. 

The information about the shell motion is again encoded in the evolution of p+. Charts 
that describe the shell motion in the exterior geometry more transparently can be con- 
structed as in section |V|. 



VIII. SUMMARY AND DISCUSSION 

In this paper we have analyzed the Hamiltonian structure of spherically symmetric Ein- 
stein gravity coupled to an infinitesimally thin null-dust shell. We formulated the theory 
under Kruskal-like boundary conditions, prescribing the evolution of the spatial hypersur- 
faces at the two spacelike infinities. We adopted smoothness conditions that made the 
variational equations distributionally well defined, and equivalent to the Einstein equations 
for this system. 

We then simplified the constraints by a Kuchaf-type canonical transformation and per- 
formed the Hamiltonian reduction. It was seen that the reduced phase space consists of 
two disconnected copies of M^, one for a right-moving shell and the other for a left-moving 
shell. We found on each component a global canonical chart in which the configuration 
variables are the Schwarzschild masses on the two sides of the shell, leaving the shell dy- 
namics indirectly encoded in the conjugate momenta. Excluding the special case of a shell 
straddling a horizon, we found a local canonical chart in which the configuration variables 
are the shell curvature radius and the interior mass, in an arbitrarily specifiable stationary 
coordinate system exterior to the shell. In particular, performing a partial reduction and 
fixing the interior mass to be a prescribed constant, we reproduced a previously known shell 
Hamiltonian in the spatially fiat gauge outside the shell. 

We also cast into canonical form the theory in which the evolution at the infinities 
is freed by introducing parametrization clocks. We found on the unreduced phase space a 
global canonical chart in which the physical degrees of freedom and the pure gauge degrees of 
freedom are completely decoupled, and we identified the pure gauge configuration variables in 
this chart as the Eddington-Finkelstein time and the curvature radius. Finally, we adapted 
the analysis to the spatial topology M^, which has just one infinity, and for which the 
spacetime inside the shell is fiat. Expectedly, the reduced phase space for this spatial 
topology turned out to consist of two disconnected copies of M?, one for an expanding shell 
and the other for a collapsing shell. 

In addition to the Kruskal spatial topology 5*^ x M and the Euclidean spatial topology R^, 
yet another spatial topology of interest would be that of the RP geon 0, RP^\{a point 



at infinity}. As the reduced phase space of the vacuum theory with the RP geon topology 
has dimension two one expects that the reduced phase space with a null shell would 
have dimension four. Indeed, this is the conclusion reached under a technically slightly 
different but qualitatively similar falloff in Ref. by first performing a Hamiltonian 
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reduction for a massive dust shell and then taking the zero rest mass limit. It does not seem 
straightforward to adapt the canonical transformation of section III to MP^ geon topology, 
however. An MP^-geon-type spacetime with a null shell can be mapped to a Kruskal-type 
spacetime with two null shells, but these two shells must be moving in opposite directions; 
our canonical transformation, on the other hand, was adapted to only one direction of the 
shell motion at a time. 

A similar issue arises if one wishes to include more than one null-dust shell. One expects 
our canonical transformations to generalize readily to the case when all the shells are moving 
in the same direction. Shells moving in different directions would, however, seem to require 
new methods. 

Several steps in our analysis relied crucially on the fact that the shell is null. This 
issue appears first in the consistency of the ADM equations of motion in section y. In a 
fixed background geometry, the equations obtained by varying the action ( p.6|) with respect 
to the shell variables must, by construction, be equivalent to the geodesic equation for 
the shell. In our dynamical equations ( |2.9|) , the pair consisting of ( |2.9^ ) and ( p.Qf ), if 
interpreted individually on each side of the shell, must therefore be equivalent to the null 
geodesic equation. The reason why the potentially ambiguous equation ( |2.91] ) turns out to be 
unambiguous is precisely that the junction is along a null hypersurface, and this hypersurface 
is geodesic in the geometries on both sides of the junction. 

Next, the fact that the shell history is null led us to the Eddington-Finkelstein time as a 
spacetime function that is sufficiently smooth to provide an acceptable momentum conjugate 
to M(r). Finally, in section |^, the null character of the shell history made it possible to 
leave the interior canonical pair (m_,p_) untouched in the canonical transformation from 
(m_|_, p+) to {r,p)- This is because the null history, when viewed from the exterior geometry, 
does not contain information about the interior mass, beyond the statement that m_ < m_|_. 

These special properties of a null shell suggest that our analysis may not be immedi- 
ately generalizable to timelike shells. For example, for a dust shell with a positive rest mass, 
already the consistency of the ADM equations of motion fails under our smoothness assump- 
tions: the variational equations corresponding to (p.8|) and (|2.9| ) only become consistent if 
the right-hand side in the counterpart of ( p.91| ) is by hand interpreted as its average over the 



two sides of the shell fS^l. However, new avenues may open if one relaxes the assumption 



that the variations of the geometry and matter be independent. Recent progress in this 
direction has been made by Hajicek and Kijowski [p^ -|5^. 

The work in this paper has been purely classical. One may, however, hope that our 
canonical charts on the reduced phase space will prove useful for quantizing the system. In 
the spatially flat gauge outside the shell, the quantization of the shell variables with fixed 
interior mass was introduced as a model for Hawking radiation with back-reaction in Ref. 
Il30| , and the same approach was applied to related black holes in Refs. Our results 

provide the tools for a similar analysis in an arbitrarily specifiable stationary gauge outside 
the shell. Whether this freedom in the gauge choice can be utilized to a physically interesting 
end remains to be seen. One may also wish to explore quantizations based on the global 
canonical charts in which the dynamics is simpler but the spacetime picture more hidden. 
This might shed light on the analogous question of quantizing in a dynamically simple but 
geometrically nontransparent canonical chart in the context of a two-dimensional dilatonic 
gravity theory coupled to scalar fields We leave these questions subject to future work. 
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APPENDIX A: HAMILTONIAN EQUATIONS OF MOTION AT THE SHELL 

In this appendix we isolate the independent information that the Hamiltonian equations 
of motion, ( p.8| ) and ( |2.9| ), contain at the shell. It will be shown in appendices ^ and ^ 
that when this information is combined to Einstein's equations away from the shell, we 
unambiguously recover the correct junction conditions for general relativity coupled to a 
null-dust shell. 

To begin, we note that equations ( ^.81 ) and all save the last one of equations ( |2.9| ) have an 
unambiguous distributional interpretation. The constraint equations ( |2.8|) contain explicit 
delta-functions in r from the matter contribution and implicit delta-functions in R" and P^. 
The right-hand sides of ( |2.9a|) and ( |2.9b|) contain at worst finite discontinuities, and the 



right-hand sides of (|2.9cD and (|2.9d|) contain at worst delta-functions; this is consistent with 
the left-hand sides of ( p.9aD -( p.9d] ), recalling that the loci of nonsmoothness in A, i?. Pa 
and Pr may evolve smoothly in t. Wherever explicit or implicit delta-functions appear, they 
are multiplied by continuous functions of r. The only potentially troublesome equation is 
therefore ( p.91| ): the right-hand side is a combination of spatial derivatives evaluated at the 
shell, but our assumptions allow these derivatives to be discontinuous. 

If / stands for any of our metric functions that may be discontinuous at the shell, we 
define 

A/:= hm [/(r + e)-/(r-e)] . (Al) 

The delta-contributions to /' and / at the shell can then be written respectively as (A/)5(r — 
r) and — r(A/)5(r — r). With this notation, the constraint equations (|2.8|) at the shell read 

AP' = -r/p/P , (A2a) 
APa = -p/A , (A2b) 

and the delta-contributions in the dynamical equations ( |2.9c| ) and ( p.9d| ) at the shell read 



- rAPA = + A^'APa , (A3a) 
A2 

iVAP' + RAN' ^ , , , , 

-rAP^ = . -I- N^APr . (A3b 

A 

The full set of equations at the shell therefore consists of ( p.9e|) , (|2.9f|) , (|A2|) , and ([A^). Of 
these, all except ( p.91| ) are manifestly well defined. 



Two of the six equations are easily seen to be redundant. First, inserting APa from ( |A2b|) 
into ( |A3aD yields an equation that is proportional to ( |2.9e| ) by the factor p/A. Equation 
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( ]A3a|) can therefore be dropped. Second, by continuity of the metric, we observe that 
R{t) = R{t, x{t)) is well defined for all t, and so is its total time derivative, given by 



R=[R + xR'i 



(A4) 



The individual terms on the right-hand side of ( |A4| ) are not continuous at the shell, but the 
left-hand side shows that the sum must be, and we obtain 



AR 



-xAR' 



(A5) 



An entirely similar reasoning leads to counterparts of (|A5|) with R replaced by any 
metric function that is continuous in r. Using ( |A5D and ( |2.9e| ), equation ( |2.5bD gives 
APa = f]R{AR')/A. This shows that the two equations in ( |A2|) are equivalent, and we 
can drop (|A2b| ) . 

To simplify ( |A3bD , we evaluate APr from ( |2.5bD and eliminate R and A using ( |A5| ) and 
its counterpart for A. Using (p.9e|), the result can be arranged to read 



= A liVaVn'] 



where the vector field v"" is defined by 



(A6) 



(A7a) 
(A7b) 



both at the shell and away from the shell. At the shell, f " coincides with the shell history 
tangent vector ( |2.10| ), by virtue of the equation of motion ( |2.9e| ). Through standard 
manipulations, ([A6|) can be brought to the form 



= A [vbV^Vaidrf] 



(A8) 



where Va is the spacetime covariant derivative. 

The information in the Hamiltonian equations of motion at the shell is therefore captured 
by the set consisting of {^^, (glQl ), (^), and (|X8|). 

In appendices and |C| we combine these four equations to the fact that away from 
the shell, equations ( p.8| ) and ( p.9| ) are equivalent to Einstein's equations and thus make 
the geometry locally Schwarzschild. As noted in subsection |IIB| , equation ( p.9e|) implies 
that the shell history is null. We therefore only need to examine two qualitatively different 
cases, according to whether or not the shell history lies on a horizon. The results, derived 
respectively in appendices ^ and y, are summarized here in the following two paragraphs: 

When the shell history does not lie on a horizon, the continuity of R across the shell 
completely determines the geometry. Equation ( [A^ ) reduces to an identity, and the com- 
bination of derivatives on the right-hand side of ( p.91| ) is continuous at the shell. Equation 
( p.91|) becomes then well defined. With p given by (|A2a|) , equation (|2.9fl) reduces to an 
identity. 

When the shell history lies on a horizon, the masses on the two sides must agree. Equa- 
tion ( |A8| ) now implies that the soldering along the junction is affine. The combination of 
derivatives on the right-hand side of (p.9f|) is then continuous at the shell, and equation 
( ^.91| ) becomes well defined. With p given by ( [A2aj ), equation ( ^.91| ) reduces to an identity. 
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APPENDIX B: EQUATIONS OF MOTION FOR A NONSTATIC SHELL 



In this appendix we verify the claims in the penultimate paragraph of appendix ^. For 
concreteness, and without loss of generality, we may assume that the shell history lies right 
of the horizon that the shell does not cross. The geometry is then as in figure ^for rj = —1, 
and its time inverse for rj = 1. 

On each side of the shell, we introduce the Eddington-Finkelstein coordinates, 

ds'^ = -FdV^ - 2r]dVdR + R^dn^ , (Bla) 
F=l- 2M/R , (Bib) 

where M is the Schwarzschild mass. To avoid cluttering the notation, we suppress indices 
that would distinguish the coordinate patches on the two sides of the shell. Wherever 
ambiguous quantities are encountered [such as in equations (|B4|) below], the equations are 
understood to hold individually on each side of the shell. 

The coordinates {t, r) of section || can be embedded in the metric ( pi|) as = V{t, r) 
and R = R{t, r), independently on each side of the shell. We obtain 

gtt = -FV^~2T]VR , (B2a) 
grr = -FV'^ - 27]V'R' , (B2b) 
g^^ = -FVV - 7] (VR' + V'r) . (B2c) 



As the surfaces of constant t are spacelike, rjV < everywhere. Expressing the ADM 
variables in terms of the metric components and using (p^ ), we find 

!^-iV^ = -^ . (B3) 
A V ^ ^ 

As in section ^ the shell history is written as r = r(t), and we write R{t) := R(t, x{t)). 
We also write, independently on each side of the shell, V(t) := V(t,t(t)), and similarly for 

R, R', V, V, and so on. We then have, independently on each side of the shell, 

R = R + iR' , (B4a) 
V = V + xV' . (B4b) 



With these preliminaries, we turn to the shell equations of motion. First, equation ( p.9e|) 



implies, with the help of (|B3|) and (P4b|) , that = 0. The shell history is therefore a 

hypersurface of constant V, independently on each side. These equations also imply that 

V/V is unambiguous and 

x = -V/V' . (B5) 



Equations ( P2b| ) and ( |B2c|) yield, after eliminating V and R with the help of (|B4aj) and ( |B5D 
the relation 
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rjV'R = -XQrr - gtr 



(B6) 



As the right-hand side of (|Bq ) is unambiguous, and as R ^ 0, (|Bq ) imphes that V is 

unambiguous. Thus, both V' and V are unambiguous. 

Consider next the constraint ( |A2a|) . As V and unambiguous, equation ( |B2b|) 

imphes AR' = —^rjV'AF. Using ( |Blb|) , the constraint ( |A2a| ) becomes 



-V'AM 



(B7) 



Consider then equation ( p. 91] ). Using 
forward calculation yields 



and the relation V = V + xV", a straight- 



V_ 

V' 



m 



As the right-hand side of ( [B8| ) is unambiguous, equations ( [B3[ ) and ( PH]) show that the 



right-hand side of ([2.91]) is unambiguous. Further, when (|B7| ) holds, it is seen that (|2.9fl) is 
identically satisfied. 

What remains is equation ([A^). In the coordinates (Bl) we have, again independently 
on the two sides of the shell, = (0, v^, 0, 0), and 



-7] 



V 

TJ ,R 



(B9) 



As = i?, is unambiguous. As (dr)^ = V (drY = V, we see that (dr)^ is unambiguous. 
As the vector field denoted on each side by d/dR is continuous at the shell and tangent to 

the shell history, is unambiguous. Therefore, the right-hand side of (p9|) is 

unambiguous, and equation (|A^) is identically satisfied. 



APPENDIX C: EQUATIONS OF MOTION FOR A STATIC SHELL 

In this appendix we verify the claims in the last paragraph of appendix 0. For concrete- 
ness, and without loss of generality, we may take rj = —1, so that the shell is moving to the 
left, and the geometry is as in figure ||. 

On each side of the shell, we introduce the Kruskal null coordinates, 

ds"^ = -G du dv + R^dn'^ , (Cla) 
G = (32MVi?) exp(-i?/2M) , (Clb) 

-uv = - l) exp(i?/2M) , (Clc) 

where M is the common value of the Schwarzschild mass. When there is a need to distinguish 
the two coordinate patches, we write the coordinates as {u±,v±), with the upper (lower) 
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sign referring to the patch on the right (left). The ranges of the coordinates are f + > and 
V- < 0, with — oo < u± < oo, and the shell history lies at the common horizon at v± = 0. 
When the index is suppressed, equations containing ambiguous terms are understood to hold 
individually on each side of the shell. 

Embedding the coordinates (t, r) in the metric (PH ) as m = u{t,r) and v = v{t,r), 
independently on each side of the shell, we obtain 

gtt = -Guv , (C2a) 
grr = -Guv' , (C2b) 
gt, = -IG {iiv' + u'v) . (C2c) 

As the surfaces of constant t are spacelike, f ' > everywhere. Expressing the ADM variables 
in terms of the metric components and using (C^f), we find 



As above, we introduce the quantities ii, ii, u', and so on, and similarly for v. The 
counterparts of equations (B4) read 

u = u + iu' , (C4a) 

{, = ^+iv' . (C4b) 

Equation ( p.9e| ) then implies, with the help of ([C^ ) and ( |C4b| ), that -0 = 0, v/v' is unam- 
biguous, and 

r = -i/v' . (C5) 



Consider next equation ( |A8|) . In the coordinates ( |C1|) we have, independently on the 
two sides of the shell, v"- = (f", 0, 0, 0), and 

[vkV'^Vaidrff = [iv'')'idr)u,uf • (C6) 

Equation ( [A 81 ) therefore reads 

{dr)u+,u^ = K")' {dr)u.,u. , (C7) 

at the junction v± = 0. Writing = {du^/du+)v'^+ , {dr)u. = {du+/du-){dr)u+, and 
du_ = {du+/du-)du_f_, we obtain d'^U-/du\ = 0. This means that at the junction t>± = we 
have 

ii+ = au_ + 13 , (C8) 

where a and (3 are constants and a > 0. As the Kruskal coordinates are afiine parameters 
along the horizons, this means that the soldering of the two geometries along the common 
horizon is affine. 

Consider next the constraint ( [A2a| ). From (|Clc|) we have R' = —2Muv'. The constraint 
( [A2a| ) therefore reads 
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-AM'^Aiuv') 



(C9) 



Equations (|C2b|) and (|C2c|) yield, after eliminating v and ii with the help of (|C4a|) and (|C5| ), 
the relation 



9tr 



(CIO) 



As the right-hand side of 



is unambiguous, v'u is unambiguous. The affine relation 
ID imphes = aul, and hence 



Hence A{uv') = and equation ( [C5| ) takes the form 

p = -4M^(3v\. . 



(Cll) 



(C12) 



As p < by assumption and > 0, we have (3 > 0. This means that the right-hand- 
side bifurcation two-sphere occurs earlier on the history than the left-hand-side bifurcation 
two- sphere, as shown in figure 

Consider finally equation (|2.9f|) . Using (|C5| ) and proceeding as with (|B8| ) , we find 



v' 
v' 



(C13) 



By ( |C11| ), the right-hand side of ( |C13| ) is unambiguous. Equations ( |C3| ) and ( |C13| ) then 
show that the right-hand side of ( p.9f ) is unambiguous. When (|C12| ) holds, it is seen that 
(p. 91]) is identically satisfied. 
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FIGURES 




FIG. 1. The Penrose diagram for a spacetime in which the shell does not straddle a horizon. 
The shell history is the dashed line passing through points pi and p2- The shell has been taken 
left-moving, which means r/ = — 1, and to lie in the future of the left-going horizon, which means 
that the right-hand-side Schwarzschild mass m+ is greater than the left-hand-side Schwarzschild 
mass m_. The diagrams for r/ = 1 and/or m_|- < m_ are obtained through inversions of space 
or time or both. The spacetime is uniquely determined by the values of m+, m_, and r]. A 
hypersurface of constant t extends from the left-hand-side to the right-hand-side ^ and the 



points A and B indicate the ends of the asymptotically null hypersurface introduced in section IV 



The dotted lines are hypersurfaces of constant null time ending respectively at A and B. Point B 
is here shown as being in the future of the shell history, but in general it could be anywhere on the 
right-hand-side X~. 
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FIG. 2. The Penrose diagram for a spacetime in which the shell straddles a horizon. The shell 
history is the line passing through points pi and p2- The spacetimes on the two sides share a 
common Schwarzschild mass m. The shell has been taken left-moving, which means ry = — 1. The 
diagram corresponding to r/ = 1 is obtained through time (or, equivalently, space) inversion. The 
spacetime is uniquely determined by the values of m and rj. The points A and B and the dotted 
null lines ending at them are as in figure |l|. 
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